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The possibility of coexistence of superconductivity (SC) and antiferromagnetic long range order 
(AFLRO) of the two-dimensional extended t — J model in the very underdoped region is studied by 
the variational Monte-Carlo (VMC) method. In addition to using previously studied wave functions, 
a recently proposed new wave function generated from the half-filled Mott insulator is used. For hole- 
doped systems, the phase boundary between AFLRO and d— wave SC for the physical parameters, 
J/t = 0.3, t' /t — —0.3 and t"/t = 0.2, is located near hole density Sc = 0.06, and there is no 
coexistence. The phase transition is first-order between these two homogeneous phases at Sc. 

PACS numbers: 74.20.-z, 74.25. Ha 



Correlation between the d-wave SC and AFLRO is 
one of the critical issues in the physics of the high- 
temperature superconductivity (HTS)P, 0. Early ex- 
perimental results showed one of the common features of 
the HTS cuprates is the existence of AFLRO at temper- 
ature lower than the Neel temperature T^r in the insu- 
lating perovskite parent compounds. When charge car- 
riers (electrons or holes) are doped into the parent com- 
pounds, AFLRO is destroyed quickly and then SC ap- 
pears. In most thermodynamic measurements, AFLRO 
does not coexist with SCQ. However, this is still a con- 
troversial issue. Recent experiments such as neutron- 
scattering and muon spin rotation show that the spin 
density wave (SDW) may compete, or coexist with SC 
under the external magnetic field [j, |5|, p . Remarkably, 
elastic neutron scattering experiments for underdoped 
YBa2Cu30^ {x = 6.5 and 6.6, = 55K and 62.7K, 
respectively) show that the commensurate AFLRO de- 
velops around room temperature with a large correla- 
tion length ~ lOOA and a small staggered magnetization 
mo ~ 0.05/.tBlllEli. These results suggest that AFLRO 
may coexist with SC but the possibility of inhomogeneous 
phases is not completely ruled out. 

For the theoretical part, the two-dimensional (2D) t- 
J model is the first model proposed|lClj to understand 
the physics of HTS. Anderson proposed the resonating- 
valence-bond (RVB) theory for the model about one and 
a half decades ago. The theory is reexamined again [llj 
recently. The authors compared the prediction of the 
RVB theory with several experimental results and found 
the theory to have successfully explained the main fea- 
tures of cuprates. This so called "plain vanilla" theory 
did not consider the issue of AFLRO, which must be 
addressed at very low doping. From analytical and nu- 



merical studies of the t — J model, it was shown that 
at half-filling, the d— wave RVB state with AFLRO is a 
good trial wave function (TWF). In this case, SC cor- 
relation is zero because of the constraint of no-double- 
occupancy. Upon doping, the carriers become mobile and 
SC revives while AFLRO is quickly suppressed. However, 
if the doping density is still small, AFLRO will survive. 
Thus SC and AFLRO coexist in the very underdoped 
regime [H H Q E [H . Exact diagonalization (ED) 
up to 26 sites show that both SC and AFLRO are en- 
hanced by the external staggered field. This result also 
implies these two orders can coexist homogeneously in 
a 2D t ~ J model^3- However, the regime of AFLRO 
predicted by these studies extend to larger doping than 
the experimental results. The robustness of the coexis- 
tence of SC and AFLRO seems to be inconsistent with 
experiments ISj. 

There are several experimental and theoretical studies 
suggesting the presence of the next- and third-nearest- 
neighbor hopping terms t' and t" in cuprates. For ex- 
ample, the topology of the large Fermi surface and the 
single-hole dispersion studied by ARPES and the 
asymmetry of phase diagrams of the electron- and hole- 
doped cuprates can be understood by introducing these 
terms. Further, these longer range hopping terms may 
be essential for the large enough Tc for the t — J— type 
models (20! l2l|. In this paper, we'd like to demonstrate 
the phase diagram constructed by VMC results of the ex- 
tended t — J model. The trial WF's for very underdoped 
systems are generalized from the siiigle-hole and slightly- 
doped WF proposed by Lee et al^^ |2^. The results 
for the hole-doped case show there is no coexistence of 
d— wave pairing and AFLRO when the next- and third- 
nearest neighbor hopping terms are introduced. And the 
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phase boundary of AFLRO is pushed to lower doping 
density. 

The Hamihonian of the extended t — J model is 



H =Ht + Hj = 



(1) 



1 



- YUj{claC],a + H.C.) + J ^ (Si • Sj + ^n^Uj) 
ij <id> 

where tij = t, t', and t" for sites i and j are nearest, 
next nearest, and the third nearest neighbors. < i,j > 
in Hj means the spin-spin interaction occurs only for 
nearest neighbors. Ci^a — (1 — ^ii,-cr)C(,cr, satisfies the no- 
double-occupancy constraint. At half-filling, the system 
is reduced to the Heisenberg Hamiltonian Hj. As carriers 
are doped into the parent compound, Ht is included in 
the Hamiltonian. 

First we exam the phase diagram of the t — J model, 
that is, t' = t" = 0. Following Ref.[il E^, three 
mean-field order parameters are introduced: the stag- 
gered magnetization = {S\) = —{Sq), where the 
lattice is divided into A and B sublattices, the uniform 
bond order parameters x — (X^o- '^L'^io')' ^^'^ d— wave 
RVB (d— RVB) one A = {cjiCii — Cjicu) if * J 
n.n. sites in the x direction and —A for the y direction. 
The Lee-Shih WF, which is the mean- field ground state 
WF is 
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where Ns is the total number of sites and Ak 
ek)/Ak and Bk - -'^''^ 



(2) 



iK"' - Ck )/Ak with 



(^k ' + Kf'^ and = + A2)i/2. Here Ak = 
|A(cosfca; — cosky). Energy dispersions for the two SDW 
bands are — ±[(ek + y^)^ + (Jm^)^]^/^ — ^ with ek = 
-2{t5 -I- '^Jx){coskx + cosky). Ok^ = akCko- + cr/3kCk+Q<T 
and 6ko- = -cr/3kCkff + "kCk-i-Qcr, where Q (7r,7r), 
< = Hi - [(^k + + m)]} and Pi = \{l + [(ek + 

/^)/(^k + ^^)W^ a,re the operators of the lower and up- 
per SDW bands, respectively, /x is the chemical potential 
which determines the number of electrons. Note that the 
summation in Eq. ^ is taken over the sublattice Brillouin 
zone (SBZ). The operator Pd enforces the constraint of 
no doubly occupied sites for cases with finite doping. 

For the half filled case, /j, = and the optimal vari- 
ational energy of this TWF obtained by tuning A and 
TOs in the VMC simulation is — 0.332J per bond which 
is within 1% of the best estimate of the ground state 
energy of the Heisenberg model. For the case of pure 
AFLRO without A, energy per bond is about 3 to 4% 
higher. Upon doping, there are two methods to modify 
the TWF: one is to use the SDW bands with a nonzero 
/i, the other is to create charge excitations from the half- 
filled ground states. For the former method, the TWF is 
optimized by tuning A, mg and /i. Note that for larger 



doping densities, AFLRO disappears (mg = 0) and the 
WF reduces to the standard d — RVB WF. For the latter 
method, the WF is the "small Fermi pocket" state | ^p): 
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k ^ Qp means the k points in the Fermi pocket Qp are 
not occupied. For example, for 4 holes in 12 x 12 lattice, 
Qp = {(7r/2, 7r/2), (7r/2, -7r/2)}. The number of holes is 
twice of the number of k-points in Qp and is identical 
to zero in Eq.©. In general, for the ground state the set 
Qp should be determined variationally. Yet as we expect, 
it agrees well with th e rig id band picture in the slightly 
doped cases as in R.ef. |23| . 

The staggered magnetization {M) = 
■ e^'^'^jS';^.) and the d-wave pair-pair cor- 



ir(E.Ai, Ar,+r), where 
0.3 and t' = t" = for the 

r>ave 
^d 



relation Pci(R) = 

Arj = CRi|(cRj+x| 

are measured for J /t 

12 X 12 lattice with periodic boundary condition, 
is the averaged value of the long-range part (| R |> 2) of 
Pd(R). The resulting {M) (fuU circles) and P^'"'' (empty 
circles) are shown in FigQ] 
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FIG. 1: (M) (full symbols) and P^" (empty symbols) for 
J/t = 0.3, t' = t" ^ (circles) and t' /t = -0.3, t" /t = 0.2 
(triangles) for hole doped 12 x 12 lattice. 

It can be seen in Fig^ that in the underdoped re- 
gion for the J/t = 0.3, t' ^ t" ^ case, AFLRO (fuU 
circles) coexists with SC (open circles) for dc < 10%, 
which is smaller than the weak-coupling mean-field re- 
sult ~ 15%|l5j, but still larger than the phase boundary 
of AFLRO determined by experiments Sc < 5%. The 
energies of | are lower than those of | for all 

doping densities in this case. This result is also consistent 
with the results reported by Himeda and OgataJ^. The 
VMC result is more realistic than the weak-coupling one. 
It may result from the rigorous no-double-occupancy 
constraint that suppresses the AFLRO faster than the 
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constraint-relaxed mean-field approximation. 

Now we examine the phase diagram for J/t — 0.3, 
t'/t = -0.3 and t" /t = 0.2. For this case, the WF 
Eq.(|2Jl is modified by replacing ^ by ^ + 4t'^cosk^cosky + 
2t'^{cos2k.j; + cos2ky) due to the second and third nearest 
neighbor hopping terms, t'^ and are variational pa- 
rameters, t'^ and t" are not necessarily equal to the bare 
values t' and t" because the constraint strongly renor- 
malizes the hopping amplitude. On the other hand, the 
effect of t' and t" on | is the choice of k-points in 
Qp, and the form of Eq.Q is not changed. 



yet P^^^ is greatly suppressed by replacing ^ by Qp to 
control the density. This replacement seems to make the 
WF decoherent for pairing. 

For 6h larger than 0.06, the RVB state (m^ = in 
I ^'ls)) optimizes the energy. P^"*^ increases and (M) 
drops to zero sharply. Unlike the t' ~ t" = case there 
is no region optimized by | ^ls) with non-zero m^. In 
conclusion, there is no coexistence of AFLRO and SC for 
the t'/t = —0.3,t"/t = 0.2 case. These parameters are 
close to the values for YBCO and BSCO comDOunds[23|. 
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FIG. 2: Energy difference per site between the two wave func- 
tions I *p> and I *LS> for t'/t = -0.3 and t" /t = 0.2 in 
12 X 12 (full circles) and 16 x 16 (open squares) lattices. 

The optimal wave functions for different densities are 
determined by minimizing the variational energies among 
^p{ms,A,{Qp})) and | ^-lsI™., A, t^, i^', /i)). The dif- 
ferences of the energies of best | ^fp) and | ^ls) for 
various hole densities are shown in FigI21 (M) (full tri- 
angles) and P^""^ (open triangles) for 12 x 12 lattice are 
shown in Fig^ 

It can be seen that level crossing occurs at Sh ^ 0.06. 

^I^p) has lower energy below the critical density. To 
show I ^'ls) and | ^fp) belong to two different types 
of WF, we calculate the overlap of them, (y^ffjpq") 
is only 0.0113(4). The almost orthogonality of the two 
wave functions implies that the ground state WF's switch 
at the critical density. Another evidence is shown by the 
correlation functions of the two wave functions shown in 
FigEl It is clear that the holes in | \l/p) repel each others 
and pairing is very small, while the behavior is opposite 
for I *Ls). 

For Sh < 0.06, I ^fp) is the ground state WF and (M) 
is a little larger than the t' = t" = case while P^^^ 
is suppressed by one order of magnitude. Thus there 
is AFLRO but no SC in this regime. The behavior is 
quite different from | ^ls) for the same doping regime 
for t' = t" = case, whose coexists with (M) . 

The possible reason is that the WF | ^fp) gains energy 
(short-range effect) from its d — RVB feature as | 
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FIG. 3: Pair-pair correlation Pd{R) (circles) and hole-hole 
correlation H{R) (triangles) of [ 'i/Ls) (full symbols) and | Vfp) 
(open symbols) for 8 holes in a 12 x 12 lattice. 

In summary, for the extended t — J model, we proposed 
a new WF | ^'p) for the underdoped regime which has 
lower variational energy than the traditional WF with 
coexisting AFLRO and SC. This WF is constructed un- 
der the framework of RVB. The new wave function has 
AFLRO but SC is largely suppressed and there is no co- 
existence of AFLRO and SC in the underdoped regime 
of the hole-doped extended t — J model. The variational 
phase diagram shows better agreement with experimen- 
tal results for the underdoped HTS cuprates. 

Note that in this study we only consider the homo- 
geneous states. Since the phase transition comes from 
the level crossing of the two classes of states at the crit- 
ical density 6c = 0.06, it is a first order phase transi- 
tion. It is quite natural to have inhomogeneity in the 
system near the critical point l^^. It may also lead to 
other more novel inhomogeneous states such as stripe 
phase 26J . Another interesting result of our study is that 
the non-coexistence of SC and AFLRO is much more ro- 
bust for systems with larger values of t'/t and t" /t such 
as YBCO and BSCOH- For LSCO where t'/t and t" /t 
are smaller, the tendency toward coexistence is larger 
and the possibility of inhomogeneous phase will become 
much more likely. 
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